disorder-induced enhancement of transport through graphene p-n junctions 
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We investigate the electron transport through a graphene p-n junction under a perpendicular 
magnetic field. By using Landauar-Buttiker formalism combining with the non-equilibrium Green 
function method, the conductance is studied for the clean and disordered samples. For the clean 
p-n junction, the conductance is quite small. In the presence of disorders, it is strongly enhanced 
and exhibits plateau structure at suitable range of disorders. Our numerical results show that the 
lowest plateau can survive for a very broad range of disorder strength, but the existence of high 
plateaus depends on system parameters and sometimes can not be formed at all. When the disorder 
is slightly outside of this disorder range, some conductance plateaus can still emerge with its value 
lower than the ideal value. These results are in excellent agreement with the recent experiment. 1 
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Due to the recent success in fabrication of graphene, a 
single-layer hexagonal lattice of carbon atoms, great at- 
tention has been attracted in the research of graphene. 2,3 
The unique band structure of graphene with a linear 
dispersion relation (E = ±7ii>|k|) near the Dirac-points 
leads to many peculiar properties^ For instance, the 
quasi-particles obey the Dirac-like equation and have the 
relativistic-like behaviors with zero rest mass. Its Hall 
plateaus assume the half-integer values g(n + l/2)e 2 /h 
with the degeneracy g = 4£& By varying the gate volt- 
age, the carrier of graphene can be tuned from electron- 
like to hole-like and vice versa. To examine the interplay 
between the electron-like and hole-like quasi-particles, a 
graphene p-n junction would be a good candidate. Many 
exciting phenomena reflecting the massless Dirac char- 
acter of carriers^* 7 - such as relativistic Klein tunneling 5 
and Veselago lensingS., were predicted for the graphene 
p-n junction. 

Very recently, the graphene junction has been realized 
experimentally^ As expected, it was found that in quan- 
tum Hall regime the two-terminal conductance exhibits 
quantized plateaus with half-integer values for the p-p 
or n-n junctions. For the disordered p-n junction, new 
plateaus emerge at e 2 /h and (3/2)e 2 /h. At about the 
same time, a theoretical analysis qualitatively explained 
the appearance of these plateaus that is due to the mix- 
ture of the electron and hole Hall edge modes in the p-n 
boundary^ After that, subsequent works have also inves- 
tigated the graphene p-n junction^ However, these the- 
ory can not accounted for the experimentally observed 
plateau that appeared at about lAe 2 /h which is lower 
than the expected value (3/2)e 2 //i. In addition, the rea- 
son that the expected plateaus at 3e 2 /h or higher values 
have not been observed remained mysterious. In view of 
this situation, it is clear that a quantitatively analysis for 
the graphene p-n junction is needed. 

In this paper, we theoretically study the electron trans- 
port through the p-n junction of disordered graphene 



under a perpendicular magnetic field B. By using the 
tight-binding model and the Landauar-Buttiker formal- 
ism combining with the non-equilibrium Green function 
method, the conductance is calculated in both clean and 
disordered samples. Numerical results show that the con- 
ductance is very weak in the clean p-n region at large 
B. With increasing of the disorder, the conductance 
is strongly enhanced and new plateaus emerge at e 2 //i, 
(3/2)e 2 //i, etc. The range of the disorder strength W 
needed for the existence of the lowest plateau e 2 /h is very 
broad, so this plateau can easily be observed. But for the 
plateaus corresponds to higher quantization values, the 
range of W can be very narrow. Sometimes these higher 
order plateaus can not be formed. Our results seem to 
suggest that the higher the plateau value, the more diffi- 
cult it is to observe experimentally. When the disorder is 
slightly off this disorder range, the conductance plateaus 
can still emerge, but its value is lower than the expected 
one. These results are in excellent agreement with exper- 
imental data. 

In the tight-binding representation, the Hamiltonian 
of the graphene p-n junction (see Fig. la) is given byJ^ 
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where a\ and a% are the the creation and annihilation 
operators at the discrete site i, and et is the on-site en- 
ergy. In the left and right leads, e; = El or Er, which 
can be controlled by the gate voltages. The disorder ex- 
ists only in the center region. The potential drop from 
the right to the left leads is assumed to be linear, i.e., 
ei = k(E R - E L )/{2M + 2) + E L + w it where M is 
the length of the center region and k = 0, 1, 2, ...2M + 1 
(see Fig. la). The on-site disorder energy Wi is uniform 
distributed in the range [-W/2, W/2) with the disorder 
strength W. The size of the center region is described by 
the width N and length M, and it has 2N(2M + 1) car- 
bon atoms. In Fig. la, it shows a system with N = 2 and 



M = 4. Here we only consider the zigzag edge graphene, 
but all results are similar for the armchair edge graphene. 
The second term in Hamiltonian (1) describes the nearest 
neighbor hopping. Due to the existence of the perpendic- 
ular magnetic field B, a phase faj is added in the hopping 
element, and faj = f: A ■ dl/fa with the vector potential 
A = (-By, 0, 0) and fa = fi/e. 

The current flowing through the graphene p-n junc- 
tion is calculated from the Landauer-Biittiker formula: 11 
/ = (2e/h)Jde T LR (e)[f L (e) - f R (e)}, where f a (e) = 
1/ exp[(e - eV a )/k B T +1] (a = L, R) is the Fermi distri- 
bution function in the left and right graphene lead. Here 
Tlr(c) = Tr[T lGT n.G a ] is the transmission coefficient, 
where the linewidth functions r Q (e) = i[£„(e) — E°(e)], 
the Green functions G r (e) = [G a (e)] t = l/[e - H cen - 
S£ — SJj], H ce „ is the Hamiltonian in the center region, 
and E£(e) is the retarded self-energy due to the coupling 
to the lead-a that can be calculated numerically^ 2 -. After 
obtaining the current /, the linear conductance is given 
by G = lim dl/dV. 

In the following numerical calculations, we use the hop- 
ping energy t rs 2.7beV as the energy unit. Since the 
hoping energy t corresponds to 10 K, we can safely set 
the temperature to zero in our calculation. The width 
N is chosen as N = 50 in all calculations. Since the 
nearest- neighbor carbon-carbon distance a = 0.142nm, 
the width is 3aN — 21.6nm for N = 50. The magnetic 
field is expressed in terms of <f> with 4> = (3\/3 / '4)a 2 B / 'fa 
and (3\/3/2)a 2 S is the magnetic flux in the honeycomb 
lattice. In the presence of the disorder, the conductance 
is averaged over up to 2000 random configurations except 
for Fig. 2b where only 400 random configurations were 
used for each data. In the experiment, the typical con- 
centration of electrons or holes is around 10 13 /cm 2 that 
corresponds to the on-site energies El, Er < O.lt. So we 
will mainly focus on the region of El and En within 0.3i. 
In this range of energy, the dispersion relation is linear 
and exhibits Dirac behaviors. 

We first study the clean graphene junction. Fig. lb and 
c show the conductance G versus the Fermi level of right 
lead En setting the magnetic field B = 0. In the n-n 
region with El, Er < 0, G is approximatively quantized 
and exhibits a series of equidistant plateaus at the half- 
integers (in the unit of 4e 2 /h) due to the transverse sub- 
bands of the lead with finite width. Because of the linear 
dispersion relation the transverse sub-bands E n of the 
confined graphene are in equidistant instead of E n ~ 
n 2 of the usual two-dimensional electron gas. While for 
Er < El, due to the fixed sub-band numbers in the left 
region, no more higher plateaus appear. On the other 
hand, in the p-n region with El < and Er > 0, there 
is no plateaus. The conductance G in the p-n region 
(En > 0) is always less than the corresponding plateau 
value in the n-n region (En < 0). Due to the occurrence 
of the Klein tunneling processes^ the conductance is quite 
large, e.g., G > e 2 /h for almost all positive Er at M = 5. 
With the increase of M, the Klein tunneling processes are 



slightly weakened and so is the conductance. 

Next, we examine the effect of the magnetic field B in 
the clean sample. With the increase of B, the equidis- 
tant sub-bands gradually evolve into the Landau lev- 
els which scales as E n ex y/n for the Dirac particle. 
The conductance plateaus in the n-n region evolve into 
the Hall plateau, and the conductance G in the p-n re- 
gion is strongly suppressed at small Er. Fig. Id, e and 
Fig. 2a show G at a high magnetic field with the mag- 
netic flux (or phase) cf> = 0.007. We see perfect Hall 
plateau in the n-n junction with equidistant in the scale 
of E R (see Fig. 2a). The plateau values are given by 
min(|^i|, \vR\)2e 2 /h where v a is the filling factors in the 
lead-a. In particular the Hall plateaus, i.e. the curve of 
G-E R in Fig.ld and le at E R < 0, do not depend on M. 
However, in the p-n region (E R > and El < 0), no 
plateaus exist. The conductance G is small and strongly 
depends on the junction length M. For the small filling 
factors vl and vr or the large junction length M, G is 
almost zero. This is because for the clean p-n junction, 
the Hall edge states for electrons and holes are well sep- 
arated in the space and do not form the mixture of the 
states leading to the very small conductance. 

In the following, we shall focus on how the conductance 
G is affected by the disorders. Fig. 2 plots G versus the 
energy El and E R with the disorder strength W = and 
W = 2. In the presence of disorders, the conductance G 
in the p-n region (or n-p region) is strongly enhanced due 
to the mixture of the electron and hole Hall edge states, 
while G in the n-n and p-p regions are slightly weakened. 
At fixed filling factors vl and v R , G is approximatively 
a constant. As vl or vr varies, a jump occurs in G with 
the borders between vl and vr regions clearly seen in 
Fig.2b. 

Now we investigate the effect of disorders on the con- 
ductance in more detail. Fig. 3 depicts the conductance vs 
E H at fixed E L = -0.1 (v L = -2) and -0.2 (v L = -6). 
When W — 0, G is small in the p-n region and G ex- 
hibits the Hall plateaus in the n-n region. With the in- 
crease of W from 0, the conductance G in the p-n region 
is strongly enhanced even for very small W . For ex- 
ample, for W = 0.02 or W = 0.05, G is greater than 
0.2e 2 /h, which is much larger than that (G < 0.001e 2 /h) 
at W = (see Fig.3a,c). When W = 0.1, the lowest 
conductance plateau with vl = — 2 and vr = 2 is well 
established with its plateau value at e 2 /h. In partic- 
ular, this plateau remains for a broad range of disorder 
strength W (from 0.1 to 3). For the higher filling factors, 
the conductance is also enhanced by the disorder, but it 
requires much large disorder to reach its ideal plateau 
value at [|^l||^'_r|/(|i / l| + \vn\)]e 2 /h. For example, for 
vl = —2 and Vr = 6 or Vl = —6 and Vr = 2, the con- 
ductance reaches the plateau of (3/2)e 2 /h when W — 2 
(see Fig. 3b, d). In the n-n region, the Hall plateau is 
not affected by the small disorders and kept their values 
at min(|/Lti|, \[j,R\)e 2 /h. If the disorder strength W is in- 
creased further, the conductance G starts to drop in both 
n-n and p-n regions. For very large W (e.g. W = 6 or 



larger) , the system enters the insulating regime and G is 
very small for all El and En. Here we wish to empha- 
size two points: (i). We have seen that the new plateau 
survives only within certain range of W — [W m i n , W max ]. 
When the disorder is slightly below W m i n or above W max , 
the conductance G still exhibits a plateau but its value 
is less than the value of ideal plateau. For example, the 
plateau of vl = —2 and vr — 6 is less than (3/2)e 2 /h 
when W = 1 and W = 3 (see Fig. 3c and 3d), (ii). For 
some high filling factor region (e.g. v L = —6 and Vr = 6), 
the conductance plateau does not emerge at all for any 
W . Because it is much more difficult to completely mix 
all states for high filling factor case, so the system goes to 
the insulating regime before the occurrence of the com- 
plete state mixture. These numerical results are in ex- 
cellent agreement with the experiment. 1 

We now focus on the conductance vs disorder strength 
for energies El and Er shown in Fig. 2a (solid dots). 
In the n-n region, the Hall edge states are very robust 
against disorders so the conductance remains quantized 
at small W (see Fig. 4b). At large disorders, the edge 
states are destroyed and the Hall conductance monoton- 
ically decreases with increasing of W . In the p-n region 
(Fig. 4a, c and d), the enhancement of conductance due to 
the states mixing at moderate disorders is clearly seen. 
For the lowest filling factors with vl = — 2 (El = —0.1) 
and vr = 2 (Er = 0.1), G reached its ideal plateau 
value e 2 /h at W = 0.09 and stayed there until W = 3 
(see Fig. 4a). We emphasize that this range of disorder 
W (from 0.09 to 3) is very broad, extends in almost two 
orders of magnitude! So the lowest plateau can easily 
be observed experimentally. For vl = — 2 and vr = 6 
(E R = 0.2) (or v R = 10 with E R = 0.25), the left side 
of the sample has an electron Hall edge state and the 
right side has three (or five) holes Hall edge states. As 
a result of the mixture of the left-side electron state and 
one of right-side hole states, the conductance G devel- 
ops a step around e 2 /h. Upon further increasing W, 
the complete mixture of the left-side electron state and 
all right-side hole states occurs at IF = 1.6 (or 1.7) 
and the conductance G reaches the ideal plateau value 
(3/2)e 2 /h (or (5/3)e 2 /h). Note that this ideal plateau 
exists only within the disorder window 1.7 < W < 2.6 
(or 1.6 < W < 2.4) that is much narrower than that of 
the lowest plateau. Our results also show that for the 



case of higher filling factors (e.g. vl — — 6 and vr = 6), 
the ideal plateau 3e 2 /h can not be reached for any disor- 
ders. When the center region becomes longer or shorter, 
the conductance G shows similar results (see Fig. 4c, d). 
For a short center region (e.g. M — 10), the conduc- 
tance reaches the ideal lowest plateau at a larger W with 
a smaller plateau width. The high conductance plateau 
at vl = —6 and vr = 6 also appears. On the other hand, 
for a longer center region (e.g. M = 40), G reaches the 
ideal lowest plateau at a smaller W with a wider plateau. 
Except for the lowest plateau, all other plateaus (includ- 
ing (3/2)e 2 /h and (5/3)e 2 /h) do not appear for M = 40. 

Finally, we study the conductance fluctuation 
rms(G) = \J (G — (G)) 2 , where (...) is the average 
over the disorder configurations with the same disorder 
strength W. Fig. 5 shows rms(G) versus W with the same 
set of parameters as in Fig. 4a and 4b. In the n-n region 
(see Fig. 5b), there is no fluctuation of the Hall edge states 
at small W. When disorder increases, the conductance 
fluctuates when the edge states are partially destroyed. 
At large disorders, rms(G) eventually goes to zero and 
enters the insulating regime. On the other hand, in the 
p-n region (see Fig. 5a), the fluctuation rms(G) is small 
for both small and large W. But rms(G) is large for in- 
termediate W and usually exhibits a double-peak struc- 
ture. In particular, rms(G) does not have the plateau 
although the conductance has a very long plateau espe- 
cially at vl — — 2 and vr = 2. 

In summary, the electron transport through a graphene 
p-n junction under the perpendicular magnetic field is 
numerically and quantitatively studied. We find the con- 
ductance is quite small for the clean p-n junction. But 
the disorder can drastically enhance the conductance 
leading to the conductance plateaus. The lowest con- 
ductance plateaus can sustain for a very broad range of 
disorder strength (about two orders of magnitude), but 
the higher plateaus are difficult to form. When the dis- 
order is slightly outside of this disorder range, some con- 
ductance plateaus in G vs Er curve can also emerge with 
plateau value smaller than the ideal value. 
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FIG. 1: (Color online) (a). The schematic diagram for a 
zigzag edge graphene p-n junction, (b), (c), (d), and (e): the 
conductance G vs Er for different center lengths M at W — 0. 
The parameters El = —0.1 for (b) and (d) and El = —0.2 
for (c) and (e), and <j> — for (b) and (c) and (j> — 0.007 for 
(d) and (e). 



FIG. 2: (Color online) The conductance G (in the unit of 
2e 2 /h) vs E L and E R with M = 20, <j> = 0.007. (a): W = 
and (b): W = 2. 



FIG. 3: The conductance G vs Er for the different disorder 
strengths W, with the parameters M = 20, <j> — 0.007. (a) 
and (b): E L = -0.1. (c) and (d): E L = -0.2. 



FIG. 4: The conductance G vs disorder strength for El and 
En fixed at different points shown in Fig. 2a with <j> — 0.007. 
The system size, M = 20 for (a) and (b), M = 10 for (c), and 
M — 40 for (d). El and En in (a) and (c) are same as in (d). 



FIG. 5: rms(G) vs Er. The parameters in (a) and (b) are 
the same as in Fig. 4a and 4b, respectively. 



1 2 



2M+1 




.4-0.3 



4-0.3-0.2-0.1 0.0 0.1 0.2 0.3 0.4 



1.0 2.0 3.0 3.6 




0.2 0.3 



1.0 

^0.8 

^0.6 



O 



0.4 - 



0.2 




(M 

CD 



2 - 



O 1 - 







1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 III 1 1 1_ 

9 1 (a)-| — ~tq ^ (b^ 



w=o 

W=0.02 

W=0.05 

W=0.1 

W=0.4 

W=1.0 



1 I \f- 

llll|llll|llll|llll P llll|ll i1l llll| 



Sp-Q 



I I I I | I I I I | I I I I | I I I 1 1 1 




-0.4-0.3-0.2-0.1 0.0 0.1 0.2 0.3 0.4-0.3-0.2-0.1 0.0 0.1 0.2 0.3 0.4 



(M 

CD 
C\J 



1.5 - 



1.0 - 



^0.5 - 



?:e: 



1.0 - 



(M 

CD 
CM 



O0.5 - 



0.0 




0.01 



0.1 



W 



10 



0.1 



w 



10 



